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ABSTRACT 

Aims. Given their photospheric origin and refractive nature, SiS molecules can provide important constraints on the relative roles of 
dust condensation and non-equilibrium processes in regulating the chemistry in circumstellar envelopes around evolved stars. 
Metliods. New SiS multi-transition (sub-)millimetre line observations of a sample of AGB stars with varying photospheric C/0- 
ratios and mass-loss rates are presented. A combination of low- and high-energy lines are important in constraining the circumstellar 
distribution of SiS molecules. A detailed radiative transfer modelling of the observed SiS line emission is performed, including the 
effect of thermal dust grains in the excitation analysis. 

Results. We find that the circumstellar fractional abundance of SiS in these environments has a strong dependence on the photospheric 
C/O-ratio as expected from chemical models. The carbon stars (C/0 > 1) have a mean fractional abundance of 3.1 x 10"*, about an 
order of magnitude higher than found for the M-type AGB stars (C/0 < 1) where the mean value is 2.7 x 10"'. These numbers are in 
reasonable agreement with photospheric LTE chemical models. SiS appears to behave similar to SiO in terms of photodissociation in 
the outer part of the circumstellar envelope. In contrast to previous results for the related molecule SiO, there is no strong correlation 
of the fractional abundance with density in the CSE, as would be the case if freeze-out onto dust grains were important. However, 
possible time-variability of the line emission in the lower J transitions and the sensitivity of the line emission to abundance gradients in 
the inner part of the CSE may mask a correlation with the density of the wind. There are indications that the SiS fractional abundance 
could be significantly higher closer to the star which, at least in the case of M-type AGB stars, would require non-equilibrium chemical 
processes. 
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1. Introduction 

Over the past decade it has become increasingly evident that 
non-equilibrium chemical processes are important in regulat- 
ing the chemistry in the circumstellar envelopes (CSEs) found 
around evolved, mass-losing, low- to intermediate mass stars 
located on the asymptotic giant branch (AGB). From equilib- 
rium chemistry calculations it is expected that the molecular 
content of the gas and dust that composes the CSEs to a large 
extent is dictated by the C/O-ratio in the photosphere of the star 
(see reviews by Glassgold 1999 and Millar 2003 and references 
therein). For M-type (O-rich) AGB stars (C/0<1) most of the 
carbon atoms will be tied up in CO leaving an oxygen dominated 
chemistry with very low abundances of other carbon bearing 
molecules such as HCN and CS. For carbon stars (C/0>1), on 
the other hand, very low abundances of oxygen bearing species 
such as, e.g., H2O and SiO are to be expected. 

However, surprisingly large amounts of oxygen bear- 
ing molecules such as, e.g., H2O, H2CO, C3O and SiO 
have been found in the CS E aro u nd the carbon star 
IRC+10216 (Keadv & Ridgwav '1993; 'Melnick et al.' '2001 
Ford et al. 2004; Hasegawa et al. 2006; Agiindez & Cernichari_ 
20061; iT^enbaum et af.ll2006t ISchoier et al.ll2006al) and, in the 
case of SiO, for a large sample of carbon stars indicating that 
such 'anomalous' chemistries might be a common phenomeno n 
in carbon-rich circumstellar envelopes (ISchoier et al.l l2006bh . 



Various suggestions to explain the observed molecular abun- 
dances have been made including non-equilibrium chemical 
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Similarly, in M-type AGB stars relatively large amounts of car- 
bon bear ing molecule s such as HCN and CS have been found 
(Olofss on et al.|[T998h . Recent interferometric line observations 
of HCN have also sorted out a long standing debate as to its ori- 
gin, clearly indicating that it is effectively formed in, or close to, 
the photosphere also in the case of M-type AGB stars jMarvell 
2005). 

It is not unreasonable to expect that also the grain type setups 
in the CSEs around mass-losing red giant star s are affected by 
these v arious chemical processes. Recently, Ho fner & AndersenI 
(I2007h have suggested that non-equilibrium processes forming 
both carbon and silicate grains in the winds of M-type AGB stars 
may help these winds to accelerate to the observed values. Such 
detailed calculations are currently not possible in hydrodynam- 
ical simulations for dust dr iven winds from M-type AGB stars 
containing pure silicates (IWoitkel 120061; iHofner & AndersenI 
l2007i) . 

Observations of SiO line emission have been shown to be 
a useful probe of the formation and evolution of dust grains in 
CSEs, as well as CSE dynamics. A major survey of SiO line 
emission from AGB stars with varying photospheric C/O-ratios 
and rn ass-l oss rates was perforin ed by iGonzalez Delgado et al.l 
(l2003h and (ISchoier et alj|2006bl) . To much surprise, detailed ra- 
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Table 1. Telescope data. 
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[GHz] 
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["] 


7 = 5-^4 


90.772 


13 


oso 


0.60 
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7 = 6^5 


108.924 


18 


oso 


0.50 


35 


y = 12 ^ 11 


217.817 
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JCMT 


0.70 


22 


y = 19 -> 18 


344.779 
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0.70 
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JCMT 


0.63 


14 


y = 20 ^ 19 


362.906 


183 


APEX 


0.70 


17 



diative transfer modelling reveal that the circumstellar fractional 
abundance of SiO has no apparent correlation with the C/O-ratio 
of the star. For the carbon stars the fractional abundance can be 
several orders of magnitude higher than predicted by thermal 
equilibrium chemistry. Interestingly, a trend of decreasing SiO 
abundance with increasing mass-loss rate of the star, interpreted 
as an effect of increased adsorption of SiO onto dust grains with 
increasing mass-loss rate, i.e., the density of the envelope, was 
found for both the M-type and carbon star samples. These claims 
have been further corroborated by recent, high qual ity, interfero- 
metric SiO observations of two M-type A GB stars (jS choier et al. 
12004.) and the carbon star IRC-H10216 dSchoier et al.. ,2006a) 
where also high spectral resolution infra-red observations of ro- 
vibrational transitions were used. 

Thus, there exist strong indications that the circumstellar 
SiO line emission carries information on the complex region, 
both chemically and dynamically, where the mass loss is initi- 
ated, and where the dust formation takes place. However, these 
conclusions still rest on somewhat loose ground, e.g., the con- 
straints on the SiO abundance distribution are poor, and the rel- 
ative importance of freeze-out onto dust grains, photodissoci- 
ation, and circumstellar chemistry is still uncertain. Important 
complementary information may be obtained by observing an- 
other Si-bearing species, e.g., SiS which has been detected in 
both oxygen- and carbon-rich envelopes (jLindqvist et alJI 19921: 
lBuiarrabalet"anil994t lOlofsson et alJll998l) . This species dif- 
fers from SiO in the sense that it can be readily formed though 
gas-phase reactions in the inner wind (Scalo & Slavskv 1980). 
Further more, it co uld possibly react differently from SiO to 
shocks (iCherchneffii2006) . dust condensation, and photodissoci- 
ation. The SiS hnes are expected to be optically thin, as opposed 
to the SiO lines which are generally optically thick also for the 
lowest mass-loss-rate objects, and thus has the potential to probe 
material even closer to the star. 

In this paper we present new (sub-)millimetre line observa- 
tions of SiS for a sample of AGB stars with varying photospheric 
C/O-ratios and mass-loss rates. The observations are supple- 
mented by a detailed non-LTE radiative transfer analysis in order 
to obtain reliable circumstellar SiS fractional abundances. The 
results are then compared to SiO abundance estimates for the 
same sample of stars and to predictions from available chemical 
models. 



2. Observations 

Multi-transition SiS line observations (see Table [1]) were per- 
formed during 1995-2006 using the Onsala 20 m telescopeQ 



' The Onsala 20 m telescope is operated by the Swedish National 
Facility for Radio Astronomy, Onsala Space observatory at Chalmers 
University of technology 



(OSO), the Swedish-ESO submilHmetre telescopeQ (SEST), the 
JCMT telescop43, and in July 2006 using the APEX telescopeQ. 

The SEST, OSO and JCMT observations were made in a dual 
beamswitch mode, where the source is alternately placed in the 
signal and the reference beam, using a beam throw of about 11' 
(SEST and OSO) or 2' (JCMT). This method produces very flat 
baselines. At the APEX 12 m telescope the observations were 
carried out using a position-switching mode, with the reference 
position located +3' in azimuth. The raw spectra are stored in 
the r* scale and converted to main-beam brightness temperature 
using r,nb - r^/'7mb- is the antenna temperature coiTected for 
atmospheric attenuation using the chopper-wheel method, and 
?7mb is the main-beam efficiency. Regular pointing checks were 
made on SiO masers (SEST and OSO) and strong CO sources 
(JCMT and APEX) and typically found to be consistent with the 
pointing model within a 3". 

The adopted beam efficiencies, together with the FWHM of 
the main beam (Omb), for all telescopes and frequencies are given 
in Table [1] The uncertainty in the absolute intensity scale is esti- 
mated to be about +20%. In Table[T]also the energy of the upper 
level involved in the particular transition (/iup) is given, ranging 
from 13 K for the 7 = 5 level up to 183 K for the 7 = 20 level, 
illustrating the potential of these multi-transition observations to 
probe a large radial range of the CSE (see Sect. |4]l. 

The data were reduced by removing a low order polynomial 
baseline and then binned (typicaly to a velocity resolution of 
about 2 km s ' ) in order to improve the signal-to-noise ratio, us- 
ing XS0 

The observed spectra are presented in Fig. [T]and velocity- 
integrated intensities are reported in Table|2] The intensity scales 
are given in main-beam brightness temperature scale (T^b). In 
addition to the new data presented here we ha ve also used SiS 
line int ensities (7=5^4 & 7=6— >5 ) reported bylBuiarrabal et al.l 
(1 19941) . lOlofsson et all (Il998h and IWoods et al.l (l2003b . In total 
1 1 carbon stars and 8 M-type AGB stars have been detected in 
SiS line emission. These stars provide our sample and are hsted 
in Tabled 

3. SiS excitation analysis 

3.1. Radiative transfer model 

The CSEs are assumed to be spherically symmetric, produced 
by a constant mass loss rate (M), and to expand at a constant ve- 
locity (ve). The SiS excitation analysis is performed using a de- 
tailed non-LTE radiative transfer code, base d on the Monte Carlo 
method and is described in more detail in ISchoier & OlofssonI 
(2001). The code has been extensively tested and benchmarked 
a gainst a wide variety of mo lecular-line radiative transfer codes 
in lvan Zadelhoff et all ( |2002|) . 

^ The SEST was located on La Silla, Chile and operated jointly by 
the Swedish National Facility for Radio Astronomy and the European 
Southern Observatory (ESO). 

^ Based on observations obtained with the James Clerk Maxwell 
Telescope, which is operated by the Joint Astronomy Centre in 
Hilo, Hawaii on behalf of the parent organizations PPARC in the 
United Kingdom, the National Research Council of Canada and The 
Netherlands Organization for Scientific Research 

* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fr Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 

' XS is a package developed by P. Bergman to reduce and analyse 
a large number of single-dish spectra. It is publically available from 
|f tp : //yggdr asi 1 . oso . Chalmers . se | 
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Fig. 1. New observations of SiS (v = 0, 7— - 1) line emission. The velocity resolution is ^2-3 km s ' in all spectra. 
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Table 2. Integrated (/obs = / T^mb dv) line intensities in Kkm s ' for the new observations of SiS (v = 0, 7— >y-l) line emission. 
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4 6^5 12 ^ 11 19 ^ 18 19 ^ 18 20 ^ 19 
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LP And 
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VCyg 

V821 Her 

CWLeo 

RWLMi 

V384 Per 

IRAS 15082-4808 

IRAS 15194-5115 
M-type 

TXCam 

RCas 

GXMon 

WXPsc 

IKTau 

IRC-10529 

IRC+40004 

IRC+50137 



0.2 : 

60 
1.2 
0.5 



0.43 



1.9 

160 



3.0 



0.80 



1.0 
0.40 



0.65 



1.5 

1.4 
4.0 
4.7 
1.3 
1.1 
1.5 



0.8 



9.2 



4.3 
1.0 



A colon (:) marks a low S/N detection. 



Table 3. Model results. 
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The SiS envelope size is determined from Eq.[3j ^^^Z^^^Z^Z^Z" 
U Cam has a detached shell that complicates the analysis of the present-day mass-loss characteristics (for details see lSchoier et al.ll2005ah . 



The circumstellar physical properties of the gas, such as 
the density, temperature, and kinematic structures, of the sam- 
ple stars is based on radiative transfer modelling of multi- 
transitional (sub-)millimetre CO line obser vations. The CO data 
used in the analysis have been presented in lSchoier & OlofssonI 



('2001'), 'Olofsson et al.' (j2002h . iGonzalez Delgado et al.1 (l2003h . 
and Schoier et al. (2006g). The kinetic temperature structure is 
obtained in a self-consistent manner from solving the energy bal- 
ance equation, where the CO line cooling is directly obtained 
from the excitation analysis. The local line width is assumed to 
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be described by a Gaussian and is made up of a micro-turbulent 
component with a Doppler width of 1.0 km s"' (e-folding radius) 
and a thermal component that is directly calculated from the de- 
rived kinetic temperature structure. 

The excitation analysis includes radiative excitation through 
the first vibrationally excited (v= 1) state for CO at 4.6 /^m and 
SiS at 13yum. Both the central star (approximated by a black- 
body) and thermal dust grains distributed in the CSE can provide 
sufficient radiation fields that start to populate the v = 1 state, 
and thereby affecting the excitation in the grou nd vibrational 
state. R elevant molecular data are summarized in ISchoier et alj 
(I2005bl) and are made publicly available through the Leiden 
Atomic and Molecular Database (LAMDAfl The collisional 
rate coefficients have been extended to include more energy 
levels as wel l as ext rapolated in temperature, as described in 
ISchoier et al](l2005bl) . An ortho-to-para ratio of 3.0 was adopted 
when weighting together collisional rate coefficients for CO 
in collisions with ortho-H2 and para-Ha. For SiS the same set 
of collis ional rate coefficient s as calculated for SiO has been 
adopted (ISchoier et al.ll2005bl) . 

The majority of sources in this study are intermediate- to 
high-mass-loss-rate objects where thermal dust emission pro- 
vides the main source of infrared photons which excite the 
V = 1 state. The addition of a dust component in the Monte 
Carlo scheme is straightforward as described in Schoi er et alj 
(I2OO2I) . The dust-temperature structure and dust-density pro- 
file are obt ained from detailed ra diative transfer modelling us- 
ing Dusty dlvezic & Elitzud[T997t ). In the modelling, where the 
SED provides the observational constraint, the dust optical depth 
specified at 10/vm, tiq, and the dust condensation temperature, 
= Tiiri), are the adjustable parameters in the;t'^-analysis. The 
effective stellar blackbody temperature, T*, is usually less well 
constrained for these intemediate- to high mass-loss-rate ob- 
j ects. The SED typ ically consists of JHKLM photometric data 
jKerschbaumll 19991 and Kerschbaum priv. com.), IRAS fluxes, 
and in some cases sub-millimetre data jGroenewegen et all 
Il993h . The total luminosity (L*) of the source is obtained from 
the period-luminosity relation of Groenewe gen & Whiteloclj 
(11996) and the distance is obtained from the SED fitting. For the 
sources where good Hipparcos parallaxes exist the correspond- 
ing distance has been adopted and the luminosity is then ob- 
tained from the SED modelling. Amorphous carbon dust grains 
with the optical constants given in lSuhl (1200 0) are adopted for the 
carbon stars while astronomical silicates are adopted for the M- 
type AGB stars. For simplicity, the dust grains are assumed to be 
of the same size (a radius, a^, of 0. 1 yum), and the same mass den- 
sity (ps = 2.0 g cm""*). The corresponding dust opacities, Ky, were 
then calculated from the optical constants and the individual 
grain properties using standard Mie theory dBohren & Hufl'manI 
Il983h . 

The parameters obtained from the CO excitation analysis and 
the dust modelling describing the physical properties of the cir- 
cumstellar envelopes are reported in Table [3] 
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Fig. 2. Upper panel - SiS fractional abundance (/o) obtained 
from an excitation analysis, as a function of a density measure 
(M/ve), for carbon stars (filled circles) and M-type (O-rich) AGB 
stars (filled triangles). Middle panel - Same as above for SiO for 
the same sample of stars (Schoier et al. 2006b; Schoier et al., in 
prep.). Lower panel - The SiS/SiO-abundance ratio. 



3.2. SiS abundance distribution 

The SiS fractional abundance distribution is assumed to be de- 
scribed by a Gaussian 



/('•) = /o exp 



(1) 



where / = «(SiS)/n(H2), i.e., the ratio of the number density 
of SiS molecules to that of H2 molecules. In circumstellar en- 
velopes such as these H is expected to be mainly in molecular 
form. In Sect.|4|we will discuss deviations from a Gaussian frac- 
tional abundance distribution. 

The best fit model is found by minimizing the total de- 
fined as 
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where / is the integrated line intensity and cr the uncertainty in 
the measured value (usually dominated by the calibration uncer- 
tainty of +20%), and the summation is done over independent 
observations. 

The size of the SiS molecular envelopes are assumed to 
be the same as those for SiO, i.e., it is ass umed to behav e 
in the s ame way to photodissoci ation (e.g., Wirsich' |1994|) . 
iGonzalez Delgad o et al.. (.2003.) and lSchoier et al.. (.2006b) found 
that a scaling law 

log re = 19.2 + 0.48 log J^j, (3) 

where M is the mass-loss rate and Ve the expansion velocity of 
the wind, provide a good fit in the case of SiO for a large sam- 
ple of M-type and carbon stars. Using this scaling law only one 
free parameter remains, the fractional abundance of SiS (/o). The 
adopted sizes using Eq.[3]are reported in Table [3] In Sect. l4.3.T] 
the validity of this approximation is tested for the two carbon 
stars CW Leo and RW LMi which are the only sources where 
high angular interferometric observations of SiS line emission 
exist. 

4. Results 

4.1. SiS abundances 

The derived abundances for the sample of 19 AGB stars where 
SiS line emission has been detected are reported in Table |3] and 
Fig. |2l It is found that the fractional abundance of SiS varies 
substantially between the carbon star (filled circles in Fig.|2]i and 
M-type samples (filled squares in Fig. |2]l. The mean fractional 
abundance of SiS in carbon stars of 3.1 x 10"^ (11 objects) is 
about an order of magnitude higher than found for the M-type 
AGB stars (on average 2.7 x 10"^; 8 objects). 

There is a weak trend that the fractional SiS abundance de- 
creases as the density of the wind (M/ve) increases. The Pearson 
correlation coefficients are r = -0.36 for the carbon stars and 
r = -0.33 for the M-type AGB stars. In comparison, for the 
same sample of sources, the derived SiO abundances have much 
stronger correlations with r - -0.60 for the carbon stars and 
r = -0.64 for the M-type AGB stars. The correlation in the case 
of SiO is further increase d for the l arger samp le of sources 
presented in Gonzalez Del gado et al.l ( l2003l) and ISchoier et al.l 
(I2006bl) and interpreted as increased importance of adsorption 
of SiO molecules onto dust grains as the density in the wind in- 
creases. 

The fractional abundances of SiS derived agree s well, typi- 
cally w ithin a factor of two, with those obtained bv Woods et al.l 
(I2003h using a simple excitation model, for the five high mass- 
loss rate objects in common. The abundances derived in the 
present analysis are usually higher than the ones reported by 
IWoods et al.l j2003h possibly a result of their adopte d excita- 
tion temperature being too high. lBuiarrabal et aTl d 19941) find SiS 
fractional abundances that on average are higher for the carbon 
stars than for the M-type AGB stars, supporting the findings in 
this paper. However, they only detect SiS emission for two M- 
type AGB stars and four car bon stars. The fractional abundances 
that iBuiarrabal et al.l (1 19941) . through a simple analysis, derive 
agree within a factor of five with the values given in Table |3]for 
the 6 sources in common. 

The fits to the SiS multi-transitional data are generally not as 
good as for CO or SiO. As examples we show the best-fit models 
for the carbon star CW Leo and the M-type AGB star IK Tau in 
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Fig. 4. Results for the best-fit model for CW Leo from the exci- 
tation analysis, a) The kinetic temperature of the gas particles as 
derived from the CO modelling is shown by the solid line. The 
dotted and dashed lines gives the excitation temperature of the 
SiS y = 5 — >4 and 7 = 20— > 19 transitions, respectively, b) The 
dotted and dashed lines gives the tangential optical depth (rt^n) 
at the line center of the SiS 7 = 5 —>4 and 7 = 20^ 19 transi- 
tions, respectively. 



Fig. |3] For CW Leo a reasonably good fit to the observations, 
both in term of the velocity-integrated line intensity (x^^i^ - 2-5) 
as well as line profile, is found. In contrast, for IK Tau both 
the line intensities (x^^^ - 9.9) and line profiles poorly fit the 
observations. In particular for the high-7 transitions, where the 
signal-to-noise is high, the models predict line profiles that are 
clearly much more flat-topped than what is actually observed. 
In Sect. |4.3| possible causes for the generally high;t'^-values will 
be discussed, including a compact, high fractional abundance re- 
gion close to the star which will naturally produce parabolic line 
profiles. 

Results from the SiS excitation analysis in the case of the 
carbon star CW Leo are shown in Fig. IH where the excitation 
temperature, Tex, and the tangential optical depth at the line cen- 
ter, Ttan, are plotted as functions of radial distance from the star 
for the 7 = 5^4 and 7 = 20 ^ 19 transitions. It is clear that both 
lines are formed under non-LTE conditions as Tex is far from the 
kinetic temperature of the gas (r^n, also shown in Fig.lUi where 
Ttan is peaking. Both lines are also optically thin (near the line 
centre) even in this high mass loss rate object. 
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Fig. 3. Best-fit models (solid lines; parameters given in TableO for the carbon star CW Leo and the M-type AGB star IK Tau over- 
layed on observed spectra (histograms). Also models with an additional compact, high SiS abundance, pre-condensation component 
are shown, (dashed line; see text for details). 



4.2. Sensitivity analysis 

In the SiS excitation analysis many parameters have been fixed 
to values obtained from modelling of other molecules (such as 
CO and SiO) or thermal dust emission. These parameters, such 
as, e.g., the mass loss rate (density structure) and temperature 
structure have some uncertainties. In Table |4] we report results 
of an extensive sensitivity analysis where the values adopted for 
several of the fixed parameters have been varied within reason- 
able limits, typically + 30% for the high mass loss rate carbon 
star CW Leo and the low mass loss rate M-type AGB star R Cas. 
The calculations were performed for a 15 m telescope (such as 
the SEST or JCMT) and the velocity-integrated line intensities 
for selected transitions from each new model were compared to 
those from the best fit model presented in Table |3] 

Varying the mass loss rate (M), i.e. the H2 density scale, by 
+ 30% results in an almost linear response in the line intensities 
indicating that the emission is more or less optically thin (see 
also Fig.m. These changes only marginally affect the line ratios 
by at most some 10%. A similar result is obtained from vary- 
ing the fractional SiS abundance (/o). The line intensities are 
significantly affected by changes in the SiS envelope size. As 
expected it is found that the lower lying (in frequency) 7 = 5 — > 4 
and 7 = 6^5 transitions, which probe material at larger radial 
distances from the star (see Fig. HJ, are most sensitive to the 
adopted SiS envelope size. The 7=19^18 and 7 = 20 — > 19 line 
emission are hardly affected at all. This means that line inten- 
sity ratios, and hence the quality of the model in fitting multi- 
transition observations, will be sensitive to the adopted envelope 
size, more so than changes in the fractional SiS abundance. 

Changing the absolute scale of the kinetic temperature struc- 
ture, obtained in a self-consistent manner from the CO analy- 
sis, by + 30% only affects the line intensities by a few percent. 
This indicates that collisional excitation of the SiS lines is not 
very effective. This is further supported by variation in the ex- 
ternal radiation fields. The excitation analysis includes the pos- 
sibility of populating also the first vibrationally-excited state of 
SiS through the absorption of 13/im photons, primarily, from 
dust emission. This can significantly change the derived line in- 



tensities of transitions within the ground vibrational state as re- 
ported in Table|4] For the high mass loss rate object CW Leo the 
thermal dust emission dominates over the direct stellar radiation 
field in the excitation. For the low mass loss rate object R Cas 
the radiation fields from the star and circumstellar dust appear 
to be about equally important. We also note here that the high-7 
lines are more sensitive to the radiation field than the low-7 lines 
affecting the line intensity ratios. 

In an analysis such as this it is very difficult to put abso- 
lute errors on the derived fractional abundances of SiS listed in 
Table |3] However, based on our sensitivity analysis a conserva- 
tive estimate would be that they are reliable to within a factor 
of two to three, given that the assumption of a single-Gaussian 
abundance profile is met. The potential of detecting chemical 
gradients in the SiS fractional abundance distribution will be fur- 
ther discussed in Sect. 14.3.21 

4.3. Wliat do the higlixii tell us? 

The resulting are, in the majority of the sources, higher than 
two, indicating relatively poor fits to the multi-transitional SiS 
data. This is in contrast to, e.g., the excitation analysis of SiO 
line emission for the same sample of sources, where x^^^ ~ 1 - 2 
in general. There are a number of potential reasons for this that 
will be addressed below. 

4.3.1 . SiS envelope size 

In the present analysis we simply assume that the SiS circum- 
stellar molecular envelope has the same size as that of SiO, i.e., 
it is assumed that these two molecules behave in the same way to 
photodissocia tion. Thi s appears to be a reasonable first approx- 
imation (e.g.. IWirsich|[r994) . However, the sensitivity tests per- 
formed in Sect. 14.21 show that the line ratios will be affected by 
the adopted envelope size and might then give poor fits to multi- 
transitional line observations if it is systematically over/under- 
estimated. High angular resolution SiS observations, capable of 
resolving the emitting regions, exist only for two sources, the 
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Fig. 5. Azimuthally averaged brightness distribution for CW 
Leo (IRC+10216), mean value around the systemic velocity 
from -24 to -28kms"' (LSR), obtained using the BIMA in- 
terferometer (data from Bieging & Tafalla 1993) as a func- 
tion of radial distance to the star in arcsec. The angular res- 
olution is 8" (FWHM). The best-fit single-dish model using 
re = 2.2 X 10"' cm and /o = 1.5 x 10"^ is indicated by the soUd 
line. A model (dash-dotted line) with re = 1.7 x 10^^ cm and 
/o =2.0 X 10"^ better reproduces the interferometer data. A still 
better fit (dashed line) is obtained by introducing a compact 
SiS component (r^ =5.0x 10'** cm and x 10"^) in addi- 

tion to the more extended envelope with re = 1.8 x 10'^ cm and 
/o = 1.7 X 10"^. All models are consistent with available single- 
dish data within Icr. 
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Fig. 6. Azimuthally averaged visibility amplitudes (real part) for 
RW LMi (CIT6), mean value around the systemic velocity from 
-2 to Okms ' (LSR), obtained using the Plateu de Bure inter- 
ferometer (data from Lindqvist et al. 2000) as a function of dis- 
tance to the phase center in kA (also shown on the upper abscissa 
is the equivalent spatial resolution in arcseconds). The best-fit 
single-dish model using re = 1.3 x lO'^ cm and /o =4.5 x 10~^ is 
indicated by the solid line. Also shown are models with a 30% 
smaller envelope size (re = 9.0x lO'^cm and /o = 8.5 x lO"''; 
dashed line) and 30% larger envelope size (re = 1.7 x lO'^cm 
and /o = 2.8 x 10"*; dash-dotted line). All models are consistent 
with available single-dish data within Icr. 
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best-fit model presented in Table|3] 



carbon stars CW Leo tee ging & Tafallalll993[) and RW LMi 
(iLindqvist et a l. 2000). In Figs. |5] and |6] we compare our best 
fit models from Table|3]with these interferometric observations. 

For CW Leo the size estimated from Eq. |3]is 2.2 x 10'* cm, 
which is clearly too large when compared with the interfero- 
metric SiS 7 = 6— » 5 BIMA observations of iBieging & TafaUal 
(Il993h as illustrated in Fig. |5] (solid line). If instead the enve- 
lope size is lowered by 23% to 1.7 x 10'* cm a significantly bet- 
ter fit to the observed brightness distribution is found (Fig. |5j 
dash-dotted line). When reducing the envelope size the frac- 
tional abundance /o was increased by 33% from 1.5 x 10"* to 
2.0 X 10 *. This new model also explains the observed multi- 
transitional single-dish data within 1 cr. 

In Fig.|6]for RW LMi the model is compared to the Plateu de 
Bure (PdB) int erferometric SiS 7 = 5 — > 4 line observations by 
ILindqvist et alj ( I2OOO1) in the (m, v)-plane in order to maximize 
both the sensitivity and the angular resolution. The best-fit model 
with re = 1 .3 X lO'* cm (from Eq.© and /o = 4.5 x 10 * (Table© 
gives a very good fit to the flux picked up by the PdB interfer- 
ometer at all baselines (Fig.|6j solid line). Also shown in Fig. |6] 
are models with a 30% smaller envelope size (re = 9.0 x 10'^ cm 
and /o = 8.5 x 10"*; dashed line) and 30% larger envelope size 
(re = 1.7 X 10'* cm and /o =2.8 x 10"*; dash-dotted line). These 
models clearly do not reproduce the observations well. All three 
models are consistent with available single-dish data within Icr. 

These two examples illustrate the usefulness of angularly 
resolved emission in further constraining the abundance dis- 
tribution by lifting the degeneracy between envelope size and 
fractional abundance that hampers multi-transitional analysis of 
single-dish data. 

Although only two sources have been studied in SiS line 
emission using interferometers, the fact that both agree within 
23% for the envelope size as estimated by Eq. [3] gives us confi- 
dence in the adopted approach, and we feel that the overall high 
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Fig. 7. Sensitivity of SiS and SiO circumstellar line emission to the inclusion of an inner compact, high fractional abundance, (pre- 
condensation) component (/c = 2 x 10"^ and r^-lx 10^^ cm) as a function of (post-condensation) fractional abundance (/o) and 
density of the wind (M/ve). The contours indicate the ratio of the total line emission for models with and without the compact 
component (see text for details). Also shown are the fractional abundances (/o) obtained from the excitation analysis for carbon 
stars (filled circles) and M-type (O-rich) AGB stars (filled triangles) with the two examples CW Leo and IK Tau high-lighted. 



X^g^ found in the SiS multi-transition analysis are generally not 
due to this adopted value. 



4.3.2. A compact, high-abundance, SiS component? 

From the fact that multi-transitional SiO excitation analy- 
ses show a clear correlation of decreasing circumstellar frac- 
tional abundance with the density of the wind (M/vg) it has 
been argued that SiO molecules effectively freeze out onto 
dust grains in the inner part of the winds from AGB stars 
jGonzalez Delgado et al.ll2003l; ISchoier et al]|2006bl) . However, 
it is difficult to put any constraints on the abundance in the 
pre-condensation region or its spatial size in the analysis of 
single-dish data as it o nly contributes ver y slightly to the total 
circumstellar emission (iGonzalez Delgado et al...2003i) . Indeed, 
high angular resolution interferometric observations are required 
in order to further constrain the abundance profiles and solid- 
ify the claim of adsorption of SiO mole cules onto dust grains 
(ISchoier et alJl2004tlSchoier et al.ll2006al) . 

Of interest is to check whether SiS line emission, as ob- 
served by single-dish telescopes in the radio regime, is more 
sensitive than SiO to such steep chemical gradients. In Fig. |7] 
the effect of including a compact region on the line intensity is 
shown as a function of post-condensation fractional abundance 
(/()) and density of the wind (M/vg). Each panel shows the ra- 
tio of line intensities of a model including also a compact, pre- 
condensation, component (with a constant fractional abundance 
of /c = 2 X 10"^ out to a cut-off radius rc = 1 x 10'^ cm) to that 



of a single Gaussian abundance distribution. A grid of models 
with varying mass loss rates and fractional abundances in the ex- 
tended component were then produced. The expansion velocity 
of the wind was take n to be 12kms~', a ty pical value found for 
AGB stars in general (iRamstedt et al.l2d06h . Thermal dust grains 
were also included in the excitation analysis with rio = 0.1. The 
star was placed at a distance of 300 pc and observations were 
simulated for SiO (7 = 2 ^ 1 and 7 = 8 ^ 7) and SiS (7 = 5 ^ 4 
and 7 = 19 — > 18) line emission using a 15 m telescope (like the 
SEST or JCMT). The spatial extents of the SiO and SiS molec- 
ular envelopes were calculated using Eq.|3]for each grid point. 

In Fig.|2]also the fractional abundances of SiS and SiO (using 
a single-Gaussian abundance distribution) found for our sample 
stars are plotted against the density of the wind (M/vg). This is 
meant as a guide to where in the parameter space the objects are 
found, and how sensitive the observed line emission is expected 
to be to a steep abundance gradient. From the location of the 
AGB stars in the panels for SiO in Fig.|7](botton row of panels) it 
is clear that any analysis based on single-dish multi-transitional 
observations is only weakly sensitive to such a compac t com- 
p onen t consistent with the findings of Gonzalez Delgad o"et al.l 
(120031) . While the line intensities may be sensitive by up to 50% 
to the inclusion of a compact component, the effect on the SiO 
7 = 8 — » 7/7 -2—^1 line ratio is less than 20% for all our sample 
sources. This means that it is generally not possible to separate 
the compact component from the extended envelope emission 
from multi-transitional analysis based on single-dish data, and 
no constraints on the compact component can be found. In addi- 
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tion, the derived post-condensation fractional abundances (using 
a single-Gaussian abundance profile) may be up to 50% too high 
in the worst case (for the majority of sources this effect is less 
than 20%). 

In contrast, for SiS it is immediately clear that the results are 
more sensitive to the inclusion of a compact component (Fig.|2l 
top row of panels). This is a combination of the lower opacities 
in the case of the SiS emission and the fact that the transitions 
readily observed by current instrumentation in the radio regime 
probe slightly warmer gas closer to the star. For the lines near 
345 GHz (see Table [B the SiS 7= 19 energy level is at 166 K, 
whereas the SiO 7 = 8 level lies at only 75 K. In particular for 
the M-type AGB stars, where the circumstellar SiS fractional 
abundances are low (~ 3 x 10"^), the inclusion of a compact high 
abundance component will drastically affect the observed line 
ratios and hence the quality of the fits to multi-transitional data 
thereby providing an explanation to the high;t'j^g^-values reported 
in Table [3] This sensitivity to chemical gradients also means 
that while there is the potential to constrain the compact, pre- 
condensation, SiS fractional abundance component (if /c is suf- 
ficiently high), the extended, post-condensation, component will 
have an added uncertainty. 

From the location of the M-type star IK Tau in Fig. [T] it 
is clear that this source may be severely affected by inclusion 
of a compact, high-abundance component. Indeed, the single- 
Gaussian abundance distribution gives a very poor fit to obser- 
vations with a reduced of 9.9 (see Table |2] and Fig.[3j solid 
line). Including a compact component with a constant fractional 
abundance of 1.1 x 10"^ out to a radius of 1.0 x 10'^ cm, 
and at the same time lowering the fractional abundance /o of the 
extended Gaussian component to 1.0x10"**, significantly im- 
proves the fit to the observations as illustrated in Fig. |3] (dashed 
line). The reduced for IK Tau using the two adjustable pa- 
rameters /c and /() is 1.8. It should also be noted that the fit 
to the various line shapes is significantly improved (better re- 
producing the parabolic shapes that are observed), in particular 
for the high-7 transitions where the signal-to-noise is high. This 
modelling of IK Tau is meant to illustrate the fact that a circum- 
stellar model including the effect of adsorption of SiS molecules 
onto dust grains can explain the observed line intensities. The ac- 
tual pre- and post-condensation fractional abundances obtained 
for IK Tau are uncertain and can not be uniquely constrained 
from the limited number of observational constraints available. 
While this procedure could easily be applied to other sources in 
our sample, thereby improving the fit to the data, more observa- 
tional constraints are needed in order to make this a meaningful 
approach. 

It has been suggested, based on observations of SiS ro- 
vibrational transitions in the infra-red, that the carbon star CW 
Leo (IRC -1-10216) has a steep SiS fractio nal abundance gradi- 
ent in its in ner wind dBovle et"aLl Il994l) . In their modelling, 
iBovle et aTl fl994) find that an initially high fractional abun- 
dance of 4.3 x 10"^, that is decreased by an order of magnitude to 
4.3 X 10"^ at a radial distance of about 8 x 10'"* cm (at an adopted 
distance of 200 pc), best reproduces the data. Observations in the 
radio regime using the Plateu de Bure and BIMA interferometers 
have also been perfo rmed (Bieging & Tafalla 1993; Lucas et al. 
ll995h . lBieging & Tafallai(ll993() find some indications of a com- 
pact SiS component from their interferometric SiS 7 = 6 — » 5 
BIMA observations. In Fig.|5](dashed line) we show that a model 
with a compact component with a constant fractional abundance 
/c of 1.7 X 10"^ out to a radius of 5 x 10'"* cm (correspond- 
ing to 8 X 10'"* cm at 200 pc) provides a somewhat better fit to 
the BIMA observations than a model without it (dash-dotted 



lin e). These va l ues ar e also fully consistent with those derived 
by ' Boyle et aTl (Il994h when considering their somewhat larger 
distance and mass-loss rate adopted for CW Leo. The location 
of CW Leo in Fig.|2]suggests that multi-transitional single-dish 
observations in the radio regime are not expected to be sensi- 
tive enough to any SiS abundance gradients. This is further il- 
lustrated in Fig. [3] where it is not possible to discriminate (the 
difference is less than the calibration uncertainty of 20%) be- 
tween a model with a compact component used to fit the BIMA 
data (dotted line) and a model using a single Gaussian abundance 
distribution (solid line) with paramete rs reported in T abled 

We note that Keadv & R idgwav (1993) and Schoie r et al.l 
(l2006al) found that explaining observed SiO line emission to- 
wards CW Leo requires a drastic decrease, by an order of mag- 
nitude, in its fractional abundance at a radial distance of about 
of 5 X lO'** cm (for a distance of 120 pc). Recent high-frequency 
observations of H'^CN 7 = 8^7 line emission towards CW Leo 
using the SMA suggest that the fractional abundance of HCN 
remains more or less constant throu ghout the inner win d and 
the extended circumstellar envelope ( ISchoier et alj|2007l) . This 
lends further support to the interpretation that both SiO and SiS 
have large gradients in their fractional abundances and that this 
is not the effect of drastically varying physical properties, such 
as density or temperature, in the envelope. Additional observa- 
tions such as these are required for AGB stars in general in order 
to obtain more reliable constraints on the SiS abundance distri- 
bution. 

4.3.3. Time variability 

It has been shown that SiS line emission can be time variable. 
ICarlstrom et al] d 19901) found that the circumstellar 7 = 5^4 
and 7 = 6^5 line emission towards CW Leo varies on a regu- 
lar basis that follows the light curve in the K-band of the star. 
It was also found that when the 7 = 5^4 line emission was 
near maximum (also maximum flux in the K-band) the 7 = 6^5 
line emission was near its minimum. The variation in line inten- 
sity between maximum and minimum is about a factor of two. 
Thus observations of these lines at different epochs, which is 
usually the case in the present analysis, will significantly de- 
grade the quality of the fit of the model to the observations. 
A plausible explanation for the time variability of certain SiS 
transitions is overlap of vibrational bands in the infra-red with 
another molecular species. This will have the effect that the 
vibrationally excited levels of SiS will be pumped by an ad- 
ditional, but selective, source modifying the excitation in the 
ground vibrational state. Recent observations of possible maser 
emi ssion in the SiS 7 = 1 1 10 7 = 14 ^ 13, 7 = 15 ^ 14 lines 
by Fonfrfa Exposito et al.l (l2006l) towards CW Leo identify ex- 
isting overlap of ro- vibrational transitions in the mid-infrared of 
SiS with other abundant molecules such as C2H2 and HCN as a 
realistic pumping mechanism. 

We note that we find no direct evidence of any maser action 
in the SiS 7 = 20 — > 19 line emission observed towards CW Leo 
using the APEX telescope. This could mean that this transition 
is not affected by the selective excitation due to the overlaps or 
that the relatively large beam of the APEX telescope, compared 
with the IRAM 30 m telescope, dilutes the emission from the 
region ve ry close to the star where th e masers are thought to 
be active. iFonfrfa Exposito et aTl (l2006h derive expansion veloc- 
ities of about lOkms ' for their observed lines indicating that 
they are mainly produced in the acceleration zone (s; 5 - 7 stel- 
lar radii) before the terminal velocity (14km s ') of the wind is 
reached. From modelling of our transitions it is clear that we 
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are not directly probing this region of the wind (see Fig. [3]). 
Currently, the effect of time variability due to maser emission 
is unknown for AGB stars in general. 

5. Circumstellar chemistry 

Stellar atmosphere models predict that the SiS fractional abun- 
dance in AGB stars is strongly dependent on the photo- 
spheric C/O-ratio (see, e.g., Millar 2003 and Cherchneff 2006). 
Typically, SiS fractional abundances under the condition of LTE 
carbon are expected to be ~ 10"^ in carbon stars whereas in M- 
type AGB-stars they are expected to be only ~ 10"^ - 10"^. From 
the results reported in Table [3] it is clear that the citcumstellar 
fractional abundance of SiS varies substantially between the car- 
bon star and M-type samples. The mean fractional abundances of 
SiS found in the present study are in reasonable agreement with 
LTE stellar atmosphere chemistries. For carbon stars the mean 
is 3. 1 X 10"^, about an order of magnitude higher than found for 
the M-type AGB stars (on average 2.7 x 10"^). 

Based on the present modelling, using a single Gaussian 
abundance profile, there is a weak trend that the SiS frac- 
tional abundance for the sample sources decreases as M/vg, 
i.e. the density, in the wind increases (Fig. |2j top panel). A 
much stronger correlation is found for the related m olecule SiO 
jGonzalez Delgado et al] 120031; [Schoier et alj |2006b) as shown 
in Fig. |2] (middle panel) for the stars also detected in SiS. This 
has been interpreted as effective adsorption of SiO molecules 
onto dust grains in the inner wind which h as been further cor- 
roborated by interfe rometric observations dSchoier et alj|2004t 
ISchoier et alj|2006al) . In addition, there appears to be no way 
of distinguishing a C-rich chemistry from that of an O-rich 
based on an estimate of the circumstellar SiO abundance alone, 
suggesting that non-equilibrium chemical processes are impor- 
tant. For SiS the direct interpretation of the derived circumstel- 
lar fractional abundances is more complicated as discussed in 
Sect. 14. 31 which could explain the weaker correlation of its frac- 
tional abundance with the density of the wind. 

A departure from LTE is to be expected by the variable 
nature of AGB stars that induces shock waves that propa- 
gates through the photosphere thereby a ffecting its chemistry. 
Models of shocked ste l lar atmospheres dWillacv & Cherchneffl 
ll998HDuari et al.lll999HCherchnefill2006h indicate that the frac- 
tional abundances of molecules, such as SiO and SiS, can 
be significantly altered by the passage of periodic shocks. 
IWillacv & Cherchnefj (Il998[) in their chemical modelling of 
CW Leo (IRCh-10216) find a thermal equilibrium value for the 
SiS fractional abundance (relative to H2) of 1.5 x 10"^ which 
is increased to 3.4 x 10' ^ due t o non-equilibrium chemical pro- 
cesses. For comparison, iMillad f2003 ) obtains an LTE SiS frac- 
tional abundance of 1.0 x 10"^. In the present analysis of multi- 
transitional SiS line emission towards CW Leo we derive a 
lower circumstellar value of 1.4 x 10"*, using a single Gaussian 
abundance component (see Sect. |4] and Table [3]l. However, as 
discussed in Sect. 14.3.21 there are observational evidence that 
the fractional abundance closer to the star may be as high as 
1.7x10"^ (see Fig. |5]) and more in line with the chemical 
model predictions. The low SiS fractional abundance in the ex- 
tended, post-condensation, part of the envelope indicates that 
SiS molecules, like SiO, are significantly incorporated into dust 
grains. 

Similarly, iDuari et akl d 19991) found a fractional SiS abun- 
dance (relative to H2) equilibrium value of 4.4 x 10"^ in their 
chemical modelling of the M-type AGB star IK Tau which 
is initially increased by an order of magnitude to 4.8 x lO"*" 



near the photosphere due to shocks but is significantly reduced 
to 8.8 X 10"'" before the dust formation zone. In comparison, 
^herchneff (2006) derived a significantly lower LTE value of 
7.1 X 10"^ for the same C/O-ratio of 0.75 and a non-equilibrium 
value near the photosphere of 1 .8 x lO"*". Also, a drastic decrease 
of the SiS fractional abundance in the chemically and dynami- 
cally active region between 1-5 stellar radii, before dust forma- 
tion sets in, is found by Cherchneff (2006) who derive a very low 
fractional abundance of » 1.0 x 10"^ for their synthetic model of 
TX Cam. In our excitation analysis we derive circumstellar SiS 
fractional abundances of 1 .0 x 10"^ for IK Tau and 4.0 x 10"^ for 
TX Cam, however, the fits to multi-transitional SiS data are poor, 
in particular for IK Tau. A significantly better fit is obtained if the 
fractional abundance close to the star is drastically increased, in 
the case of IK Tau up to 1 . 1 x 10"^ assuming a size of this com- 
ponent of 1.0 X 10'^ cm (see Sect. 14.3. 2l and Fig.[3]l. Such a high 
fractional abundance in the inner wind would favour the impor- 
tance of non-equilibrium processes. Presently, it is not clear if 
the very low fractional abundance (1.0 x 10"*^) found in the ex- 
tended envelope around IK Tau is a result of non-equilibrium 
chemistry in the inner wind or a result of freeze-out onto dust 
grains. Interferometric observations are clearly required in or- 
der to make further progress and better constrain the chemical 
models. Also the estimated ch emical reaction ra tes involving Si 
might have to be re-examined (ICherchnefIll2006h . 

6. Conclusions 

A systematic and extensive excitation analysis of multi- 
transition (sub-) millimetre SiS line emission has been per- 
formed for a large sample of 19 AGB stars with a range of photo- 
spheric C/O-ratios and mass-loss rates. By constraining the cir- 
cumstellar SiS abundance distribution we reach the following 
conclusions: 

- In order to derive reliable fractional abundances of SiS, exci- 
tation through IR ro-vibrational transitions, needs to be taken 
into account. This will in most cases (for our sample sources) 
mean that the effect of dust grains needs to be included in the 
excitation analysis. The present work contains a proper treat- 
ment of this. 

- The adopted SiS envelope size will also affect the quality of 
the abundance estimates and it has been tested against avail- 
able interferometric SiS line observations for the two carbon 
stars CW Leo (IRC-H10216) and RW LMi (CIT6). It is found 
that the empirical formula derived from SiO observations is 
adequate to use also for SiS, i.e., these two molecules appear 
to behave similarly to photodissociation. 

- The mean fractional abundance of SiS in carbon stars of 
3. 1 X 10~* is about an order of magnitude higher than found 
for the M-type AGB stars (on average 2.7 x 10"^). Taken at 
face value, these circumstellar fractional abundances are in 
reasonable agreement with photospheric LTE chemical mod- 
els. However, there are indications that the SiS fractional 
abundance could be significantly higher closer to the star 
which, at least in the case of M-type AGB stars, would re- 
quire non-equilibrium chemical processes. 

- No clear trend that the SiS fractional abundance decreases as 
the density (M/vg) of the wind increases is found, in con- 
trast to the results for SiO for the same sample of stars. 
This could mean that SiS molecules are less likely to ad- 
sorb onto dust grains than SiO molecules. However, possible 
time-variability of the line emission in the lower J transi- 
tions and the sensitivity of the line emission to abundance 
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gradients in the inner part of the CSE may mask a correla- 
tion with the density of the wind. Evidence supporting the 
presence of steep circumstellar abundance gradients are pre- 
sented for the carbon star CW Leo and the M-type AGB star 
IK Tau, which can naturally be explained by effective freeze- 
out of SiS molecules onto dust grains at distances of ~ 5 stel- 
lar radii. 

While it is clear that SiS line emission is a useful tool in 
investigating the conditions occuring in the chemically and dy- 
namically active region, where dust grains are formed and the 
stellar wind is accelerated, many uncertainties remain. In partic- 
ular, the constraints on the detailed SiS abundance distribution 
are poor In order to make significant progress, interferometric 
observations of both low- and high-/ SiS line emission in or- 
der to more firmly establish the relative importance of freeze-out 
onto dust grains, photodissociation, and circumstellar chemistry 
are required. We note that for a typical stellar distance of 300 pc, 
next generation interferometers such as ALMA will be required 
in order to fully resolve the region over which SiS could effec- 
tively freeze-out (angular size of » 0.2 - 0.4"). 
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